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Abstract 



The fine structure interval of P states in hydrogenlike systems can be de- 
termined theoretically with high precision, because the energy levels of P 
states are only slightly influenced by the structure of the nucleus. Therefore 
a measurement of the fine structure may serve as an excellent test of QED 
in bound systems or alternatively as a means of determining the fine struc- 
ture constant a with very high precision. In this paper an improved analytic 
calculation of higher-order binding corrections to the one-loop self energy of 
3P and 4P states in hydrogen-like systems with low nuclear charge number 
Z is presented. The method of calculation has been described earlier Q and 
is applied here to the excited P states. Because of the more complicated na- 
ture of the wave functions and the bound state poles corresponding to decay 
of the excited states, the calculations are more complex. Comparison of the 
analytic results to the extrapolated numerical data for high Z ions |p serves 
as an independent test of the analytic evaluation. New theoretical values for 
the Lamb shift of the P states and for the fine structure splittings are given. 
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I. INTRODUCTION 



Evaluations of the radiative corrections in higher order for bound states are an involved 
task because of the appearance of a multitude of terms, and because of the difficulties 
associated with bound state formalism. In this paper, we present an improved calculation 
of higher order corrections to the one-loop self energy of an electron in an excited 3P or 4P 
state. 

For the contribution 6Ese of the one-loop radiative correction to the Lamb shift of a 
bound electron, we have the following non-analytic expansion in powers of Z times the fine 
structure constant a 

SEs^^^J-^F. (1) 



where 



F = ^4,1 \n{Za)-' + ^4,0 + (Za) A^^o + 



{Zaf [^6,2 hi\Za)-^ + ^,1 hi{Za)-^ + A^^ + {Za) GsejJ • (2) 

The remainder function Gse,7 is of order 1 and is comprised of the terms A-j^ and higher 
coefficients. Corrections ^4^1, A^^q and Aq^2 vanish for P states. The terms A^^ (see, e.g., 
0]) and ^6,1 are known analytically. The term A^^ contains the Bethe logarithm which 
has been evaluated to 12 significant figures Results have not been obtained for ^6,0 

coefficients. In this paper, we present an evaluation of the ^5,0 coefficients for the 3Pi/2, 
3P3/2, 4Pi/2 and 4P3/2 states. The results lead to improved values for the Lamb shift of the 
respective states and to a new theoretical value for the fine structure splitting. We give an 
explicit formula for the fine structure of the 2P, 3P and 4P states as a function of the fine 
structure constant a, which can be used to obtain a value of a from experimental data. 

In this paper, we briefiy compare some of the methods that have been developed for the 
treatment of the one-loop problem. We give a brief account and illustrate the usefulness of 
the e-method for analytic evaluations. We then describe the evaluation of the high- 
energy part to the self energy, with a focus on details of the integration procedure. We then 
proceed to the low-energy part. Results of the calculation are given, specific contributions 
are discussed in detail. 



II. VARIOUS METHODS OF TREATMENT OF THE ONE-LOOP SELF 

ENERGY 

Using units in which h = c = 1 and = 4710;, we can write the integral corresponding 
to the one-loop self energy of an electron bound in a Coulomb field, 

6EsE= lim -le^ [ — [ Dl'Hk^ M) (^'h^ , „ ^ kttYI^) 

-{ij\6m{Mm, (3) 

where D'^^f,{k'^, M) is the Pauli-Villars regularized photon propagator (in Feynman gauge, we 
have D'lJik^, M) = -g^^{l/P - l/{k^ - M^)). The term 6m{M) in Eq. (|) is the one-loop 
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mass counter term as a function of M, Sm{M) = a (3/(4 vr)) m (ln(M^/m^) + l/2). 'i/j = ip^ 'y'^ 
denotes the Dirac adjoint. It is straightforward to derive Eq. (|]) with the Feynman rules of 
QED. By rescahng all variables to the electron mass scale 



muo 



we have 
(5£'sE 



m k', 



p ^ mp 



(4) 



duj' 

-(VS|5m(M')|V)- 



1 



A;'2 - M'2 



{^P\Y 



1 



^' - 1 - 7OI/' 



77m I 



(5) 



We will use variables rescaled to the electron mass in this paper, and suppress the prime of 
the rescaled variables in the sequel. Note that in our system of units, we have (e. g.) for the 
Bohr radius of the atom aeohr = l/{Za). By contrast, in atomic units, which are used for 
example in |Q , we would have the Bohr radius of length unity. 

The analytic properties of the propagators determine the location of the poles in the 
integrand in Eq. (^ as indicated in Fig. |T]. The original Feynman prescription calls for 
integrating the photon energy along the contour Cp. For the actual evaluation of the Lamb 
shift, however, a different contour of integration is used by most authors. Taking advantage of 
the analytic properties of the integrand and of Jordan's lemma, one can change the Feynman 
contour in the complex plane without changing the result of the calculation. We compare 
here the contour used in Bethe's original derivation of the Lamb shift, the contour used by 
Mohr in |9|-[T2|, and the contour used in Pachucki's e method, which is used in this paper. 

Mohr's and Pachucki's methods both depend on a division of the calculation into a low 
and a high-energy part. Mohr's method relies on the contour Cm in Fig. His low-energy 
part is determined by the part of the contour Cm where Re(u;) < e. The residues of the 
poles of the photon propagator only contribute to low-energy part in this case. It can be 
shown that the low-energy part is given by the formula 



AEl = lim 

5^0+ 



a 



a 



vr At{^ Jk<E„ 

1 



k 



Hjj — En ~ k — i6 



(6) 



(cf. Eq. (3.8) in is the Dirac Hamiltonian) . This contribution contains terms of lower 

order in {Za) than {Za^. The spurious lower order terms cancel when the low and the high- 
energy parts are added in that method. The high-energy part is obtained by Wick-rotating 
the Feynman contour for uo integration along the line with Re(ci;) = E^. In the non-relativistic 
limit, expression (H) corresponds (up to the term a/n En) to what would be expected to be 
the self energy of the electron in terms of traditional second order perturbation theory due 
to transverse modes of the electromagnetic field. 



AEf = Re(A^z.)--i?n 

TT 



£ P / d'kT ^ [ 

tf Jk<K .^^Am^J 



d\ ^ , ^ , (7) 



A=l,2 



En 



En' - k 
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where 



E eUk)ei(k) = 5^^- 



A=l,2 



(8) 



K in Eq. (0) is an appropriate energy cutoff to make the expression finite (for a derivation 
of Eq. (|^ cf. |[13|, Eq. (7-112) to Eq. (7-115) ibid., where in the relativistic case one has 
to substitute cx for 1/zV). In Mohr's method, K corresponds to En- Bethe's derivation of 
the Lamb shift, which gave the right scahng of the effect of the self energy and correctly 
identified v44 i, but did not include the contribution to the Lamb shift of order a/T\:{ZaY 
due to the anomalous magnetic moment, comprised the expression given in Eq. ([7p, but 
with a major modification. Bethe subtracted from (|^ the contribution that would modify 
the energy (or mass) of a free electron due to its self interaction. This contribution would 
give a contribution to the rest mass of any electron, and thus would be unobservable. In our 
terminology, Bethe's non-relativistic (NR) expression would be written 



5E 



2,NR 



-P 



2a 
Sir Jo 



K=m 



dkk{(j)\ 



p 



m 



Hs - (E^ -k) k 



^1 
m 



(9) 



where Hs is the Schrodinger Hamiltonian, cj) is the non-relativistic wave function, and the 
subtracted term —1/k in the integrand corresponds to the portion of mass renormalization 
attributable to the low-energy part. By taking the principal value (P), we identify the real 
part of Eq. (|^) as the energy shift, whereas the imaginary part corresponds to the decay 
width of the state Using the subtraction, Bethe disposed of the spurious lower-order 
terms and obtained a finite expression. 

In Pachucki's method (see Fig. an expression similar to Eq. (|^) is obtained in the 
non-relativistic limit for the low-energy part, but with an upper cutoff epsilon for the pho- 
ton energy. This cutoff epsilon separates the low and the high-energy parts. In the dipole 
approximation exp(zk ■ r) ^ 1, one obtains the expression 



5El = -— / dkk{ 
i-n Jo 



P 



p 



m Hs — {E^ — k) m 



10) 



(10) 



for the low-energy part in leading order. The renormalization term —1/k is gone, and the 
upper cutoff has been changed from K = m to K = e. The justification for leaving out the 
renormalization term is intimately linked to the special series expansion prescription used 
by Pachucki. 

Pachucki's method relies on the fact that the low-energy part and the high-energy part 
may formally be regarded as functions of the fine structure constant a and the cutoff pa- 
rameter e. Their sum, however, the self energy of the electron 6E, 



5E{a) = EL{a, e) + Enia, e). 



(11) 



does not depend on epsilon, provided the high and the low-energy parts are expanded first 
in a, and then in epsilon (the order of expansion plays a crucial role in that case). 

Another important point in Pachucki's method is that the spurious lower order terms 
which were present in Mohr's calculation vanish in the limit e — > 0, so we do not need to 
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take them into account. For example, in Mohr's calculation, the first spurious term {oi/ti) En 
originated from a trivial integration J^f" dka/n = (ce/Tf) En- In Pachucki's method, we would 
change the upper limit of integration to e and calculate lim^^Q Jq dk a / tt = 0. That means 
by choosing the e prescription, we not only make the expression for the low-energy part 
separately finite, but also dispose of the spurious lower order terms. That is the principal 
reason why Pachucki's method is well suited for the analytic calculation of higher order 
corrections to the one-loop self energy. 

The choice of epsilon remains arbitrary to a certain extent (it has to be because we ana- 
lytically expand in e and thus require arbitrariness). However, we must put some restraints 
on the magnitude of e. In the high-energy part, we expand the propagator of the bound 
electron in powers of the binding field V. We initially assume a fixed value for e which pre- 
vents infrared problems, but since eventually e 0, this expansion is regarded as a formal 
expansion that is not necessarily convergent. However, e may not be arbitrarily large. If we 
let e > 2m, we enclose poles not only from the photon propagator, but also from the negative 
spectrum of the Dirac-Coulomb propagator, which would significantly alter our expression 
for the low-energy part. It is also required that in the entire domain of the low-energy part, 
an expansion of the expression 



exp(ik • r) 



in the matrix element 



in powers of k • r corresponds to an expansion in powers of Za (this requirement justifies 
the so-called dipole approximation, in which we replace exp(ik ■ r) by unity to obtain the 
lowest order contribution to the self energy). The order of magnitude of r is l/{Za) in 
natural units. Thus we require e < (Za). The dominant contribution is then determined by 
the region in which the photon energy k = u — 0{{Zaf'), so that k • r = 0{Za). In this 
paper we consider the relativistic corrections up to relative order {Zaf. This corresponds to 
expanding exp(ik-r) up to (k-r)^. Our restrictions on the magnitude of e do not compromise 
the validity of analytic expansion in the parameter e. 



III. THE HIGH-ENERGY PART 



The high-energy part of the radiative correction is given by 



Eh — —ie^m 



dio r d^k 



Ch 2tt J (27r)3 
-{■ilj\6m{M)\il:), 



1 



A;2 A;2 - 



^- - 1 - 701/ 



7m I "0) 



(12) 



where we have used the Feynman gauge for the photon propagator {Dn^{k) — —gnu/k^) and 
the Pauli-Villars regularization prescription 



1 



1 



1 



k'^ + iS k'^ + iS k'^-M^ + iS' 



(13) 
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Note that we may leave out ie prescription when integrating along Ch since we take the 
difference of the integrand infinitesimally above and below the real axis on Ch- Along the 
positive real axis, the integrand has branch cuts due to the photon and electron propagators 
as depicted in Fig. |I[ The expression given in Eq. (|12D for is infrared divergent. In the 
evaluation, we start by calculating the matrix element 



1 



m — 'y^V 



(14) 



up to the order of (Za)^. As outlined in this can be achieved by first expanding the 
electron propagator in powers of the binding Coulomb field V. This leads to a 3-vertex, a 
double- vertex, a single- vertex and a zero- vertex part. The expansion can be diagrammatically 
represented as in Fig. ^. The resulting expressions are subsequently expanded in powers of the 
spatial electron momenta p*. This procedure is feasible for P states because up to order {Za)^, 
all of the resulting matrix elements converge. After performing the algebra of the Dirac 
matrices, the resulting matrix elements on the P state are evaluated by symbolic procedures 
written in the computer algebra system Mathematica For the evaluation, we first 
expand the wave function (given by the exact solution to the Dirac-Coulomb equation) 
in powers of (Za), then we apply operators in coordinate space representation and finally 
integrate the resulting expressions with the help of a set of rules that apply to standard 
integrals. The integrands which are to be evaluated for the matrix elements have lengths of 
up to 2,000 terms. 

We use a parametric representation of the mass counter term to allow for local cancella- 
tion of the divergences. It can be shown that 



Sm{M) = —ie^m [ — / 
JCh 27r J 



(27r) 



1 



k2 



UJ 



2 - k2 - M2 



2(uj + l] 



2-k2-2cJ 



(15) 



is a suitable parametric representation of the mass counter term along the contour Ch- The 
portion of mass renormalization along the contour Cl vanishes in the limit e — >^ 0, and we 
have 



6m{M) = —ie^m j — j 



a 



3 m 



duj f d^k 
Cf ^tt j (27r)3 

In (M^) + 1/2 



1 



A;2 + ie k^-M^ + ie 



2(uj+l] 



a;2-k2-2a; 



(16) 



Therefore, by (locally) subtracting the expression 

2{uj + l) 



5m 



before the final du rf'^k [1 / (cj' 



Cj2-k2-2cj 

•2 k2) — 1/(^^2 — k2 — M2)]-integration in Eq. (|T2D , we can 
subtract the divergences associated with mass renormalization. 

Note that 5mi contains the matrix element which is state dependent. Using the 

virial theorem for the Dirac-Coulomb equation ((ck ■ p) 
where is the dimensionless Dirac energy of the state %p. 
We give here the result for the renormalized matrix element 



{V)), we have {^i)) = E^, 
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-Pren = -P " StTIi 



up to (Za)^ in terms of k and uj. We have for the 3Pi/2 state, 

Pren(3Pl/2) = {Zaf 



8k2 + 4k^ + 6k2a;-3k^a; + 12a;2- lOk^o;^ 

3n 



27 k^ + 2u;-u;^ 



-6u;^ + 6k2u;^ + 6u;^-3a;^ 

+ {Za)^ - 128k^ + 48k^ + 24k^ + 96k^cj + 2k6cj-9k^a; 

+256 k^ - 156 k^ - 90 k^ - 168 k^ uj^ + 46 k^ u;^ 

+36 k^ Lu^ - 16 cj^ + 240 k^ cu^ + 126 k^ lj^ + 72 lj^ 

-98 k^ uj^ - 54 k^ - 132 - 78 k^ uj^ + 50 cu^ + 36 k^ cu^ 



324 k2 + 2u;-a;^ 



+18a;^-9u;'^ 

+ {Za)^ - 1044992 k^ + 516224 k^ + 319032 k^° + 32340 k^^ 

+587776 k^uj + 1716736 to + 26320 k^ u - 65310 k^^ to - 8085 k^^ 
-358400 k^ + 4218368 k^ - 3461600 k*^ - 1372476 k*^ 
-177870 k^° + 1469440 k^ - 9185344 k^ uj^ + 1027600 k^ uj^ 
+357630 k*^ uj^ + 48510 k^° + 1075200 - 9354240 k^ uj^ 
+ 10195136 k^ co^ + 2348976 k^ co^ + 404250 k^ co^ - 3978240 
+ 15638560 k^ cj^ - 3365600 k^cu^ - 777420 k^ - 121275 k^a;^ 
+7869120 uj^ - 12272960 k^ - 2002392 k^ - 485100 k*^ 
-8571360 uj^ + 3543120 k^ cj^ + 839580 k^ + 161700 k^ 
+5023200 uj^ + 852600 k^ + 323400 k^ - 1231440 
-450870 k'^ oj^ - 121275 k^ - 145740 - 113190 k^ a;^° 

+96390 uj^^ + 48510 k^ cj" + 16170 a;^^ - 8085 a;^^ 
612360 (k^ + 2 a; -a;'" 



for the 3P3/2 state, 

P'ren(3P3/2) = {Zaf 

-6a;^ + 6k^a;^ + 6a;^ 
+{ZaY 



2 , .2 



-8k^ + 4k* + 6k^cu -3k*tu + 12cu^ - lOk^tu 

3 



3cx;^ 



27 k^ + 2u; 



■ a; 



128k^ + 32k^ + 8k^ + 32k^u;-26k^u;-3k^u; 

+192 k^ uj'^ - 68 k^ u;^ - 30 k^ uj'^ - 56 k^ uj^ + 42 k^ u;^ 
+12 k^ - 112 a;^ + 42 k^ a;^ + 24 a;^ - 6 k^ uj^ - 18 k^ uj^ 

+36 o;^ - 2Q>k^uj^ - lOa;^ + I2k^ uj'^ + 6u;^ - 3a;^ 



324 



+{Zaf - 4179968 + 516608 k^ + 224616 k^° + 12180 k^^ 

+2351 104 k^ a; + 243712 k^ a; - 1022336 k^u;- 1031 10 k^° a; 
-3045 k^2 u - 1433600 k^ u'^ + 8960000 k^ + 667456 k^ o;^ 
-539868 k^ - 66990 k^° u"^ + 4157440 k^ - 2845696 k^ 
+1703632 k^ + 370230 k^ + 18270 k^° + 4300800 o;^ 
-17297280 k^ oj^ - 3980704 k^ a;^ + 135408 k^ + 152250 k^ o;^ 
-10752000 o;^ + 10612000 k^a;^ - 619136 k^o;^ - 449820 k^a;^ 
-45675 k^ + 15408960 uj^ - 99680 k^ o;^ + 586824 k^ 
-182700 k^ uj^ - 10647840 o;^ + 216720 k^ uj'^ + 159180 k^ 
+60900 k^ o;^ + 2896320 - 543480 k^ + 121800 k^ 
-278880 cu^ + 65730 k^ cu^ - 45675 k^ to^ + 136500 cu^° 
-42630 k^u;^^ - 42210 u;^^ + 18270 k^a;^^ + 6090 a;^^ 

N 7-, 



-3045 (v^^ I 2449440 (|k^ + 2 a; - u?^ ' 



for the 4Pi/2 state, 



-8k2 + 4k^ + 6k2a;-3k^u; + 12a;2- lOk^a;^ 
-6a;^ + 6k2(^3 + 6a;^-3a;^ 48 (^k^ + 2 - a;^^ 
+(Zq;)^ [ - 2944 k^ + 1072 k^ + 520 k^ + 2080 k^ a; - 10 k^ 



00 



5760 k^ 



3412k'^cj^ - IQSOk'^'w^ 



3640 k^a;^ 



780 k^ Lu^ - 560 + 5040 k^ cu^ + 2730 k^ lo^ 



-195 k^ 
+ 1050 k^cj' 

+ 1560 ~ 2070 k^ - 1170 k^ - 2700 - 1690 k^ uj^ 
+1030 a;^ + 780 k^ a;^ + 390 a;^ - 195 tu' 



15360 (^^ + 200-00^ 



(19) 



+ (Zq;)^ - 41518080 k^ + 20285056 k^ + 12019560 k^° 

+1140300 k^^ + 21790720 k^ a; + 67517184 k^u- 1295056 k^ u 
-2875250 k^° a; - 285075 k^^u; - 13189120 k^a;^ + 164921344 k^o;^ 
-138538016 k^u;^ - 50829380 k^ u;^ - 6271650 k^" cj^ 
+53975040 k^ - 357889728 k^ oo^ + 50610672 k^ oj^ 
+14327250 k^ oo^ + 1710450 k^^ + 39567360 cu^ - 357683200 k^ u;^ 
+404953920 k^ + 81291280 k^ + 14253750 k^ oo^ 
-147302400 + 601024480 k^ uj^ - 153435296 k^ cj^ 
-28556500 k^ - 4276125 k^ + 294387520 - 485016000 k^ o;^ 
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-59093160 k^o;^ - 17104500 k^a;^ - 325403680 a; V 160218800 k^cu^ 
+28458500 k^ uJ^ + 5701500 k^ o;^ + 198315040 uj^ + 17532200 k^ 
+11403000 k^a;^ - 56099120 o^'^ - 14180250 \<? - 4276125 k^ a;'' 
-920500 - 3991050 k2a;^° + 2826250 a;^^ + 1710450 k^a;^^ 

51609600 (k2 + 2 



+570150^;^^- 285075 a;^3 
and for the 4P3/2 state 



Pren(4P3/2) = {Zaf 



-8k^ + 4k^ + 6k^a;-3k*a; + 12a;^- lOk^o;' 



-6u;'^ + 6k^u;'^ + 6u;^ -3u;^ 



48 k^ + 2u;-u;^ 



- 2944 k^ + 752 k^ + 200 k^ + 800 k^ a; - 570 k^ u 



-75 k® + 4480 k^ uo"^ - 1652 k^ tu^ - 750 k^ cj^ - 1400 k^ uj^ 



970 k^ u'-" + 300 k^ - 2480 + 240 k^ + 1050 k^ uj" 



+600 u'' - 230 k^ - 450 k^ 00^ + 660 
-170 c;^ + 300k2cu^ + 150 tu^ - 75 cu^ 

5l 



650k2cu^ 



15360 (^k^ + 2a; 
+{Zaf 



41518080 k^ + 5273472 k^ + 2372328 k^° + 132300 k^^ 

+21790720 k^ a; + 3438848 k^ a; - 10087504 k^ a; - 1069810 k^° uj 
-33075 k^^ uj - 13189120 k^ o;^ + 88747008 k^ uj"^ + 4454624 k^ a;^ 
-5916484 k^ uj"^ - 727650 k^° a;^ + 37847040 k^ 0;=^ - 29007552 k^ uj^ 
+18082288 k^ uj'^ + 3870930 k^ uj'^ + 198450 k^° uj'^ + 39567360 u;^ 
-163681280 k^ - 34772416 k^ + 2115344 k^ 
+1653750 k^u;^ - 98918400 u;^ + 100285920 k^o;^ - 8639904 k^cu^ 



-4785620 k^a;^ - 496125 k^ 



uo 



143286080 o;^ - 1740480 k^tu^ 



+5429592 k^ o;^ - 1984500 k^ - 99145760 uj^ + 3382960 k^ cu^ 
+1829380 k^ cu^ + 661500 k^ (J + 26784800 uj^ - 5400920 k^ a;^ 



+1323000 k^ u"" - 2737840 tu^ + 563430 k^ u"" - 496125 k^ u" 



+1400140 o;'^ - 463050 k^ a;'^ - 408310 u;'' + 198450 k^ uo 
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+66150 o;^^ -33075 o;^^ 



51609600 k^ + 2 u; 



(20) 



Having calculated P, we finally integrate along Ch to obtain the result for Eii^ 



Eh — —ie m / -— 



duo f d'^k 



Ch '2'JT J (27r 



1 



1 



(21) 



(22) 



Note that as Pren = P — ^''tt-i, both P as well as the local mass renormalization term Snii are 
properly integrated with the regularized photon propagator. 
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The final integrations with respect to the photon momenta are done in a different way 
for the terms in Pren which require regularization and those which do not. The terms which 
require regularization are integrated covariantly by Feynman parameter techniques and a 
subsequent Wick rotation. These terms are not integrated along Ch, but rather along C'fj. 
They are not infrared divergent, so we may put e = for these terms. Those terms which 
do not require regularization are integrated in an essentially non-covariant way. The d^k 
integration is carried out first, then we proceed to the duj integration. 

The integration procedure for the terms which require regularization is as follows. We 
isolate in (P — 5mi) those terms which would be ultraviolet divergent if integrated with the 
unregularized photon propagator. We denote these terms by P^^^. We then evaluate 



5E, 



div 



du r d^k 



c„ 27r J (27r)3 - 



pdiv 
ren 



All terms in P^^^ require regularization. 

By simple power counting arguments, it can be shown that P^^^ is exactly the sum of 
those terms of order fc" where n > —2 for large k. Therefore, the terms contributing to Pf^^ 
can easily be isolated. The terms in Pdiv; can obviously be written as the sum of terms of 
the form 

where Pi is a polynomial in |k| and u and deg(pj) > n^ — 2. The entirely covariant integration 
procedure for the divergent terms will be outlined here. The terms in Eq. ( PHj ) need to be 
multiplied by the factors l/(co'^ — k^) and l/(co'^ — k^ — M^) from the regularized photon 
propagator. We use Feynman parameters in the form 

1 , nx"-^-^ 

dx -r— — r (24) 



AB^ Jo [A{l-x) + Bx]'^ 

to join the denominators. Identifying Ai = —lo^ + k^, A^ = —uj^ + k^ + and B = 
k."^ + 2uj — uj"^ , we have for the contribution from the unrenormalized photon propagator 
l/(u;2-k2), 

Ai {l-x)+Bx = -P + Di (25) 

where 

k = {u, k) = (cij — X, k) and Di = x^, (26) 
and for the contribution from the renormalization part of the photon propagator 

A{1- x) + Bx = -P + P>2 (27) 

where 

D2 = x^ + M (1- x). (28) 
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The energy shift due to the divergent terms is then proportional to 



where D represents either of the terms Di or D2 and it is understood that the contribution 
from these two terms must be subtracted to obtain the final result. Performing the Wick 
rotation 



we have 



UJ — > lUJ, 



-e = uj' + k' = kj + D 



where ke is the Euclidean 4- vector kg = {u!,k). We then obtain the energy shift due to the 
divergent terms proportional to 



k'^^+D 



The (straightforward) angular part of these integrals can be done by parametrizing the 
Euclidean 4-space a.s u = kg cos 7, k^ = ke sin 7 sin 6' cos0, k"^ = kg sin 7 sin^ sincj), 
ks = ke sin 7 cos 6'. We then have u = ke cos 7, |k| = ke sin 7. For the average over the 
4-dimensional angle, we utilize the formulae 

2 _ 1 f d^l^^ 4 _ 1 f d^l^^ 6 5 



^ cos 7 = -, / — cos 7 = -, / — cos 7 = — . 

27r2 ' 4' J 27r2 ' 8 J 27r^ ' 64 

The remaining radial part of the integrals can be evaluated with the help of the formulae 

1 1 
"^^^ (fc2 + D)P ^ 2D/3-2 (/5 - 1) (/3 - 2) ^^^^^ 

and 

It is easy to prove that the dependence on the temporary upper cutoff A disappears when the 
contribution due to the unregularized and the regularization part of the photon propagator 
are subtracted {Di and D2). 



As the final step, we integrate over the parameter x introduced in Eq. and subse- 
quently investigate the resulting expression in the limit M —>■ 00. The respective expressions 
vanish as M ^ 00 for all states considered in this paper. This fact is intimately linked 
to our using the entirely covariant Feynman parameter approach for the integration of the 
divergent terms. If we had used a non-covariant scheme of integration, as in then we 
would have had to take into account finite correction terms to obtain the correct result for 
the Lamb shift. 
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Note that the (divergent) spurious terms of order which are present in all of the 

matrix elements Pren, vanish after we have performed the d'^k integration in the way outlined 
above, which includes final expansion in the e parameter. 

The terms in which are finite when integrated with the unrenormalized photon 
propagator, do not need regularization. It is easy to see that can be written as the sum 
of terms of the form 



P 



fin 



E 



9j(|k|, 



[k2 + 2UJ 



UJ 



21 



(29) 



where qj is a polynomial in |k| and u whose degree is less than 2nj — 2 to insure ultraviolet 
convergence. For the 4Pi/2 state, e. g., P^^ is given by 



P 



fin 



4 k^ + 3 k^ u; + 6 cj^ 



24 k2 + 2a; 



UJ 



1472 k^ + 536 k^ + 1040 k^ - 5 k^ + 2880 k^ uj"^ 

-1706 k^ - 1820 k^ uj^ + 525 k^ uj^ - 280 uj^ + 2520 k^ uj^ 
+780 uj^ - 1035 k^ uj^ - 1350 o;^ + 515 uj'^ 

7680(k2 + 2cj-cj' 

+ {Zaf - 20759040 k^ + 10142528 k^ + 6009780 k^° 

+ 10895360 k^ + 33758592 k^ - 647528 k^ - 1437625 k^° uj 
-6594560 k^ cj^ + 82460672 k^ cu^ - 69269008 k^ cj^ - 25414690 k^ uj"^ 
+26987520 k^cj^ - 178944864 k^cj^ + 25305336 k^tu^ 
+7163625 k^cj^ + 19783680 cj^ - 178841600 k^cu^ + 202476960 k^cj^ 
+40645640 k^ uj^ - 73651200 uj^ + 300512240 k^ - 76717648 k^ uj^ 
-14278250 k^ uj^ + 147193760 - 242508000 k^ uj^ - 29546580 k^ uj^ 
-162701840 cu^ + 80109400 k^ uj'^ + 14229250 k^ cu^ + 99157520 uj^ 
+8766100 k^ uj^ - 28049560 uj^ - 7090125 k^ cu^ - 460250 uj^^ 



+ 1413125CU 
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25804800 (k2 + 2 cu-cj^ 



We have to calculate 



5Eh = —ie^m 



duj r d^k 1 



c„ 27r J (27r)3 tu^ - k2 



Sfin 



The quantity Fh defined by 



_a {ZaY 
t^H — —m — I'H 



(30) 



TC W 



(31) 



in Eq. (|^) may be expressed as 



12 



J Ch 



(32) 



where 



1 

2Tii 



d\k\ 



^ren(^,|k|). 



k2 



One can dispose of the factors in the numerator of the integrand using the following 
procedure. First write as 



Y -2uo 



SO y = k^ + 2u; — corresponds to the denominator in Eq. (p9|). The resulting expression 
is subsequently expanded. The powers of Y cancel, and the result does not carry powers of 
k in the numerator any more. 

The integrand in can be written as the sum of terms of the form 



1 



(33) 



k2 (k2 + 2UJ 

with suitable coefficients aj. The (i|k| integration can then be carried out using the formula 



1 
2tH 



d\k\ 



1 



1 



-1 



Tjidh 



2sgn(Im(cu)) (cu^ + fi) 
where sgn is the sign function defined as 



+ 



(n 



n~l 



(34) 



sgn(x) 



1 a; > 
-1 X <Q 



We identify 



Q = 2uj-uj^ 



(35) 



(36) 



in order to carry out the integration in Eq. (^). The formula Eq. ( |5^ deserves some com- 
ments. We are integrating along the real axis. Because u and Q both have an infinitesimal 
imaginary part all along the contour of integration Ch, the positions of the poles of the 
integrand are well defined. 

The branch cuts of the function J-'{uj) can be readily identified. Due to the term 
sgn(Im(ti;)) cij in the first term on the right hand side in (0), there is a branch cut along 
the positive real axis. This branch cut is caused by the photon propagator. Due to the term 
a/Q in the second term on the right hand side in (0) there is also a branch cut along the 
line where the expression Q = 2uj — u'^ assumes negative real values. This branch cut ex- 
tends from uj = 2 along the positive real axis to u; = oo. It is caused by the Dirac-Coulomb 
propagator (see Fig. |I|). 

It can be explicitly checked that the function JF(co') satisfies the equation 
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(37) 



We can divide the contour Ch in an upper contour which extends from e+iO+ to oo+iO^ 
and a lower contour which extends from oo — iO^ to e — iO"*". We then have due to Eq. 



poo pe 

= / did J^{uj + iO+) + / duJ^iuj -iO+) 

/oo roc 
dujJ^{u + iO+) + J du J^{uj + iO+y 

= [ du T(uj) + c.c. (38) 

We restrict ourselves therefore to the upper contour and understand that the complete 
result is the sum of the integral along the upper contour plus its complex conjugate. 
We then perform a change of variable to proceed to the final du integration. Defining 

u = -j^ -• (39) 



and 



we have 



Note that 



W(«) = ^^(a;), (40) 



^_ff = n7(Za)^ / duU{u). (41) 

Ju(Ch) 



Re(2 uj — + i 

u = " = for u e C^, Re(cj) e [e, 2[, (42) 

/Re(2 Lu — cu^) — iuj 



whereas 



jRe(uj^ -2uj) - io 

u = \ = for uo e C^, Re(cj) G [2, oo), (43) 

/Re(cj2 -2cj) + u 



where the argument of the square root is written in such a way as to represent a positive real 
quantity in the two above equations. So u{uj = 0) = 1, u{uj = 2) = —1, ^(00 + 20"*") = 0, and 
u G [—1, 0) for u G Cj^, Re(a;) G [2, 00). So for u G C|^, Re(a;) G [e, 2), u has a nonvanishing 
imaginary part, whereas for u G C^, Re (a;) G [2,cxd), m is a real quantity. The mapping 
u; M is one-on-one for u G C'^. 

Note that if we had chosen the lower contour Cj^, then ^(oo — i 0~^) = —00. In that case, 
the above substitution would not have had the desired property -u ^ for — 00. 
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On C^, oj and can be expressed as functions of u according to 



OJ — — - 



2 u ' 
2 {-u) 

i. e. \/Vt extends along the negative imaginary axis for Re(a;) > 2, w e [—1,0]. The result 
for U{u) (4Pi/2 state) is 



U{u) = {Zaf 



1 



1 



192 96(-1 + m) 

47 23 

+ 



+ 



23 



. 20480 10240 (-1 + m)^ 10240 (-1 + m)^ 
79 13 13 113 



+ 



+ 



40960 (-1 + -u)^ ' 2048 (l + ii)^ 2048 (1 + -u)^ ' 24576 (1 +-^)^ 



141737 



39173 



39173 



41287680 20643840 (-1 + uy 10321920 (-1 + uy 

2228617 3485527 4685519 



206438400 (-1 + m)^ 206438400 (-1 + -u)^ 330301440 (-1 + -u)^ 



499 



343 m 



+ 



1267 



46080 (-1 + u) 614400 65536 (1 + uf 

3801 11765 2715 

+ 



+ 



65536 (1 + m)^ 131072 (1 + m)^ 32768 (1 +m)^ 
69651 28021 89779 



1310720 (1 + uf 1310720 (1 + uf 15728640 (1 + uf 



The next step in the calculation is the u integration, 

i(e+iO+) 

where 



Fh = n^liZaf / duU{u) + c.c. 

Ju{e+i 0+) 



,V3/2 



It is useful to define 



e = 1 - u{{e + i 0+) = -i V2^ + e + ^ + 0(e)^/^ 

2v 2 



(44) 



(45) 



(46) 



(47) 



Note that the e-prescription calls for carrying out the li-integration from u{e + iO"*") to 
and subsequently expanding the result in powers of e up to e°. For those terms which are 
finite when integrated from 1 to we may carry out this integration without regarding the 
dependence on e. For instance. 
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/ du- - = + 0(v^), (48) 

so in the limit e — > 0, the e dependent term vanishes. Terms of the form 1/(1 — m)", however, 
introduce a divergence in 1/e. We have 

/■O 1 gl-n _ 1 

/ dw- - = — . (49) 

A(E+io+) (\-uY n-\ ^ ' 

We then add to the result of this integration the complex conjugate and subsequently expand 
in powers of e. This procedure is illustrated with some examples. For the terms proportional 
to 1/(1 — n)^ we have 

du- — 3— - = 1-1, (50) 



«(e+iO+) (1 — U)^ e ' 



SO 



0,1 11 

u{e+jo+) (1 — uy e e* 

= 2 



1 



-zy2i + e + 0(e)3/2 i + eO(e)3/2 
' i + 0(e)i/2) - f- ' +i + 0(e)V2' 



2- l + 0(e)^/2 = l + 0(e)i/2. (51) 



For terms proportional to 1/(1 — m) we have 

Ju(e+iO+) (1 — U) 



du ^ . = In (e) , (52) 



so 



du h c.c. = In (e) + In (e) 

u{e+iO+) 1 — U 



In (-i + 0(e)) + In (z + 0(e; 
2 2 



^ +2^ + 21n(V2i) +0(e) 



= ln2 + lne + 0(e)^/2 (53) 
Note that for constant terms it results 

/ du const. + c.c. = -2 const. + 0(e)^/^ (54) 

0+) 



Using the results Eq. (|51| ) and Eq. (|54|) , it is easy to show that the spurious terms of order 
{Zaf in the expression (^41) vanish after the final w-integration. 
The final result for the high-energy-part (4Pi/2-state) is 
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1/2) 



24409 499 „ , , 23 

(In 2 + Ine) 

86400 720^ ' 90 e 



(55) 



We give here the complete results for the high-energy-parts of the other states treated in 
this paper 



1/2) 



+ {Zaf 



6191 268 



24300 405 



(In 2 + Ine) 



20 

8U 



(56) 



FHm/2) = - + {Zaf 



67903 148,, , , 20 

In2 + lne) 

194400 405^ ' 81e 



(57) 



FHm,^) = ^ + {Zaf 



31399 137,, ^ , , 23 

(In2 + lne) 

86400 360^ ' 90e 



(58) 



IV. THE LOW ENERGY PART 

The low-energy part of the energy shift originates from low-energy virtual photons. The 
energy of the photons is comparable in magnitude to the binding energy of the electron 
(order {Z af). Therefore it is impossible to expand the electron propagator in powers of the 
binding field. We have to treat the binding field non-perturbatively. An expansion in powers 
of [Z a) is accomplished by considering the spatial momenta of the virtual photon and the 
electron momenta as expansion parameters. 

Choosing the Coulomb gauge for the photon propagator, one finds that only the spatial 
elements of this propagator contribute [|T|. The tu-integration along Cl is performed first, 
which leads to the following expression for E^, 



Er 



-e'P 



k\<e (27r)3 2|k| 



\a e 



j k ■ r 



Hd — {E^ — 1^. 



j k ■ r I 



[00 



|k|). (59) 



if/) denotes the Dirac-Coulomb-Hamiltonian H£, = a-p+j3m + V, 6^ is the transverse delta 
function, and a* refers to the Dirac a-matrices. The principal value of the above integral is 
the real quantity corresponding to the energy shift in one-loop order. The imaginary part 
of the Cl integration, which leads to the decay-width of the state, has been dropped in Eq. 
(^). In the matrix element 



1 



Hd — {E^ — 0) 

we introduce a unitary Foldy-Wouthuysen transformation V 

1 



(60) 



P'^ = {Ui)\{Ua'e'''-'U+) 



U{Hd-{E^-u:))U+ 



{Ua^e-'^-''U+)\Uij). 



(61) 



The lower components of the Foldy-Wouthuysen transformed Dirac wave function ip vanish 
up to {Zaf, so that we may approximate \U^p) by 
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\U^) = 10) + with = 0, 



(62) 



where |0) is the non-relativistic (Schrodinger-Pauh) wave function, and \6(j)) is the relativistic 
correction. 

We define an operator acting on the spinors as even if it does not mix upper and lower 
components of spinors, and we call the odd operator odd if it mixes upper and lower com- 
ponents. The Foldy-Wouthuysen Hamiltonian consists of even operators only. For the upper 
left 2x2 submatrix of this Hamiltonian, we find the result 



HF^ = UiHD-iE^-uj))U+ = m + Hs + 6H, (63) 
where Hs refers to the Schrodinger-Hamiltonian, and 6H is is the relativistic correction, 

SH = + (5 r + / o , q- ■ L (64) 

It is interesting to note the reason why we can ignore the lower 2x2 submatrix of the 
FW Hamiltonian in our scheme of calculation. The lower 2x2 submatrix contains the terms 
—m—{E^—uj) ~ —2 m as the dominating {Z a)° contribution. Therefore the integral vanishes 
in the limit e ^ 0. This can be checked by considering the integral in Eq. (^9|), inserting 
a spectral resolution for the Dirac Coulomb propagator and performing the integral over 
d'^k after suitable angular averaging. The upper 2x2 submatrix has no {Z a)^ contribution, 
because terms proportional to m and cancel. This submatrix contributes to the Lamb 
shift. Now we turn to the calculation of the Foldy-Wouthuysen transform of the operators 
a* exp (k • r). The expression U a* exp (ik ■ r) f/+ is to be calculated. Assuming that = |k| 
is of the order 0((Za)^), we may expand the expression U e^^ '^ in powers of (Za). 
The result of the calculation is 

UaV^ 'U+ = a'(l+i(k ■r)--(k ■ r?) ^—p' (a ■ p) (65) 

V 2 /2 m"' 



+7°^(l + .(k.r)-i(k.r)^) 

-7° Ap^p' - a ^ (r X S)^ 
2m-^ 2m^ r-^ 

^ 7° (k ■ r) (k X S)' - — 7° (k X S) 



2 m 2 m 

In the limit e ^ the odd operators in the above expression do not contribute to the self 
energy in (Za)^ relative order, because up to (Za)^ relative order, these operators only 
join the upper components of the wave function with the lower components of the Dirac 
Coulomb Hamiltonian. This contribution vanishes, as described. So one can neglect the odd 
operators. By using symmetry arguments, it can be shown easily that the last term in the 
above expression (proportional to k x S) does not contribute to the Lamb shift in (Za)"^ 
relative order for e ^ 0, either. 

Because we can ignore odd operators, and because the lower components of the Foldy- 
Wouthuysen transformed wave function vanish, we keep only the upper left 2x2 submatrix 
of Eq. ((|65D), and we write U e'^^'^U^ as 
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[/aV^ "[/+ ~ — fl + z(k -r) - -(k • r) 



(66) 



1 a 



2 m? 



2 m? 

+-^(k.r) (kxa)\ 
2 m 



^ r X cr 



This can be rewritten as 



Ua' (?^'U 



m 



(67) 



where 5?/* is of order (Za)^. It is understood that the term ^e^^'^ is also expanded up to 
the order (Za)^. Denoting the Schrodinger energy hj E {E = —{ZaYm/n^ for nP states) 
and the first relativistic correction to E by 5E, we can thus write the matrix element P^^ as 



+ 



m 



Hs-{E-Lu) + 6H- 5E 
We now define the dimensionless quantity 



p) 
— e 
m 



ik ■ r 



2 



10 + 50). (68) 



(69) 



Up to {Zay, we can write the matrix element P as the sum of the contributions (([70|), ([TTD, 
([7^), (0), (|7^), ([75|)). The leading contribution (the "non- relativistic dipole") is given by 



''°-^3!^<*' g.-(g-u,) 



The other contributions to P are |T[ 
• the non-relativistic quadrupole, 



1 



nq 



17- (0b' e 
6m 



ik ■ r 



Hs-{E- u\ 



■p e 



~ik ■ r I 



(70) 



(71) 



the corrections to the current a* from the Foldy-Wouthuysen transformation, 

1 



Psy = 6^''^ {<P\6y' 



Hs-{E- u 



■ p' e 



ik ■ r I 



(72) 



the contribution due to the relativistic Hamiltonian, 



(73) 



the contribution due to the relativistic correction to the energy, 

1 



SE 



(E-uj] 



■P 



(74) 
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and due to the relativistic correction to the wave function, 



Almost all of the above contributions are calculated using a coordinate space representation 
of the Schrodinger Coulomb propagator given in and [^]. Formulae given in [|I^ and 
T6| prove useful for the summation over the intermediate quantum numbers. For the non- 



relativistic quadrupole contribution, however, we use the momentum space representation 
due to Schwinger ||19| 



G(p, p', n) = AnrnX" 



le 



2 sin TIT 



/ ^PP T 

Ji dp 



X- 



X2(p - p')2 + (p2 + X2)(p'2 + X2) (1 - p)V(4p) 



2 ' 



(76) 



where X = y/—2mQ, t = ma/ X. For non-integer r, one may replace the complex integration 
around the origin (in the positive sense) by a much simpler integral 



le 



2 sin vrr 




dpp ^ g{p) 



dpp ^ 9{p)- 



(77) 



We then perform the calculation of Pnq in momentum space using formulae given in the 



paper by Gavrila and Costescu [|l^]. The momentum space wave functions of P states are 
given in |T9|. Calculations for Pnq become increasingly complex. It should be noted that in 
the case of the 4P wave function, one has to deal with intermediate expressions of up to 
20,000 terms. We do not describe these calculations in any further detail. Evaluations for 
this part of the calculation were done on IBM RISC/6000 systems with the help of the 



computer algebra system Mathematica 



The final u integration is done by change of variable 



UJ 



t, where t 



(78) 



(t = corresponds to = oo, and t = 1 corresponds to = 0). As an example for a 
P matrix element, we give here the result for a contribution to Psh{4:P 1/2), caused by the 
Russell-Saunders coupling term in the Foldy-Wouthuysen transformed Dirac Hamiltonian: 



Pl.s(4Pi/2) = W 
= {Za) 



Hs - {E^ - 00) [Am'^r^ 



a 



cr L 



16384 F366(t) (-1 + 2 1) (1 + 2 1) (-1 + 4 1) (1 + 4 1) (-5 + 9 1 



2\2 



(675(-l + t)2(l + t) 
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8192 <^23{t) (-1 + 2 1) (1 + 2 1) (-1 + 4 1) (1 + 4 1) (-5 + 9 1 



2\2 
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675(-l + t)'(l + t) 



14 



-{Zaf 

x(-5 + 9t 
-{Zaf 
-{Zaf 



2\2 



16384 *63(^) (-1 + 2 (1 + 2 (-1 + 4 (1 + 4 1) 

675(-l+t)2(l + t)^^^ 
16 F4(i) t^{-5 + 9 (5 + 9 i^) ^ 
40967t^(l + 2t) (l + 4t) (-5 + 9t^^^ 



45(-l + i)5(l + i)^ 



X 



X 



- 630 + 2988 1 - 22677 + 322539 - 1308681 + 5929071 

+20884413 + 53752275 - 124047785 + 233683811 
-362131987 + 489361437 - 552966475 + 52677470U^^ 
-433556447 + 298335889 - 171601869 + 84114837 
-33425148 + 10588500 - 2811264 + 536448 
-64512^^^ + 9216 



675 (-2 + i) {-l + tf{l + tfA{-2,^2t) (-7 + 4^) (-5 + 4i) (-3 + 4^) 



+ {Zaf 



54385222500+ 13984771500 1+ 15573367761525 



+ 137704964619600 



1186676073750825 - 891514989328950 



-18753057125404875 + 257731557911828220 - 644138132939443685 
-162189845176905698 + 5703957333855251853 



11 



-40461889136476779376 1 
-470465622621377811686 1^^ 
-285989274333023499852 1^^ 



12 



179879463449806219483 1 
f 722355119079050313441 
- 2287627466701513790137 



+8721622433626707698710 1" - 17878002921800012898021 1'^ + 
+22584212519300184263992 - 11732105237171234433515 
-22234628025915862014322^21 + 73751337137966874205423 
-118367926111467660583708 + 122509973591672089703233 
-68708753123836196116438^2^ - 27553239135523122581409^^6 
+126820257213436783378720^27 - 189954960243232498928039^^8 
+199312481720191257908238 P - 164901212087123003552709 
+112286946900261699896044 - 64112830789916243409479 P 
+30834230557497197343734 - 12455975426801658077884 
+4211137005166537183048 - 1 183248309588246704096 
+270125031800068986496 - 48193573673016712704 
+6583393931443034112 i^'^ - 691086520295792640 

+41007381492105216 + 1314475331420160 i^^^ 
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where 



425250 (-2 + tf (-1 + (1 + tf^ (-3 + 2tf 

x(-l + 2t) (-7 + AtY (-5 + 4t)2 (-3 + 4:tf (-1 + 4t) 



8192 r (-1 + 2 1) (-1 + 4 1) (-5 + 9 ty (3 + 6 1 + 3 r + 8 



xln(2/(l + t)) 



675(-l + t)i^(l + t) 



16 G4(t) (-5 + 9 1^) (-1125 + 2400 1 + 145390 



-19200 7 r - 148320 - 2026412 T + 418560 7 T 
+1027200 + 7316050 - 1920000 7 - 1382400 

-7444143 1*^ + 2211840 7 1*^ - 19200 ln(2/(l + t)) + 418560 ln(2/(l + t)) 
-1920000 ln(2/(l + t)) + 2211840 ln(2/(l + t))) 

10125 (-l + t)^(l + t)* 

+ {Zaf 256 (-5 + 9 (8 + 45 1 - 305 f - 175 + 283 

-265 + 3421 + 395 f + 433 1^) ln[(2 t)/(l + t)] ^ 

675(-l + t)^2(l + t)^l, (79) 



F^it) = 2F1 1, -4t, 1 - 4t, {{t - l)/(t + 1))^ 



(80) 



G^{t) = 2F1 [1, -4t, 1 - 4t, (t - l)/(t + 1)] , 



(81) 



(82) 



^"^'^"S (3-4t + ^)^ 



2fc 



(83) 



*63(t) 



E 

fc=6 



(3 - 4t + k) 



^(A; + 6), 



(84) 



where \1/ denotes the logarithmic derivative of the F function. The quadratic singularity in 
the result for Pl s(4Pi/2) in Eq. (|79|) at t = 3/4 (given by the (-3 + 4t)^-term in the de- 
nominator of the purely rational function) corresponds to the decay into the 3D state. One 
can check explicitly that the insertion of corresponding intermediate states in the spectral 
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decomposition of the propagators necessitates the existence of quadratic singularities in the 
P matrix elements, and that the quadratic singularities occur only in the matrix elements 
with two propagators. The quadratic singularities are a consequence of the perturbative 
treatment of 5H in the propagator l/{Hs + — {E^ + SE — u)) (expansion in 6H is not 
allowed in the vicinity of a pole of the resolvent G{E) = l/{Hs — E)). The integration pro- 
cedure for the quadratic singularities is as follows: first we isolate and calculate analytically 
the integral of the term that gives rise to the singularity (as a function of t), then we take 
the difference of the edge terms at t = 1 and t = 0. This procedure takes back the effect of 
the perturbative treatment and assigns the correct value to the t integral. As the final step, 
we subtract the term that gave rise to the quadratic singularity and proceed with the rest 
of the terms in the usual way described in 

The integration procedure deserves some further comments. We also encounter in the 
matrix elements singularities of linear type att = l/4, t = l/2 and t = 3/4, which also cor- 
respond to the decay of the excited state. The residue taken at these singularities yields the 
decay width of the respective states. In order to obtain the principal value of the t-integral, 
one has to symmetrize the integrand around all the singularities. This is also accomplished 
by symbolic procedures written in the computer algebra language Mathematica. 

The results of the calculations have been checked in many ways. An important cross- 
check is the cancellation of e-divergent terms in the sum of the high and low-energy parts. By 
considering the expansion of the propagators in powers of the logarithmic singularities 
of all contributions can be calculated individually and agree with the results obtained from 
complete evaluation. 

The contributions to the low-energy part Fl in {Zay relative order are given in Table 
I and Table ||. Summing all contributions, we obtain the following complete results for the 
contribution of the low energy parts F^. 



Fl{3P^/^) = --In ko{3P) + {Zaf 



(85) 



i^L(3P3/2) = --lnA;o(3P) + (Z«)2 



0.69356(1) + — + — In- ' 



lie 405 (Za)2 



(86) 



FL(4Pi/2) = -77lnfco(4P) + (^«) 



23 499 , e 
0.99780 1 + + — lu' 



90 e 720 {Zay 



(87) 



and 



Pl(4P3/2) 



■In fco(4P) + (Za) 



23 137 e 
0.73057 1 + + In- 

^ ' 90 e 360 {Zaf 



V. RESULTS AND EVALUATION OF THE LAMB SHIFT 



Summing the contributions from the high and low-energy parts, we obtain the following 
results for the scaled F function defined in Eq. (0): 
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i^(3Pi/2) = -\ - ^lnA;o(3P) + [Zaf 

6 



-1.14768(1) + ^ln[(Za)-^] 



(89) 



F(3P3/2) = ^-^lnA;o(3P) + (Z«)^ 



-0.59756(1) + l^ln[(Za)-2] 



(90) 



F(4Pi/2) = -\ - ^lnfco(3P) + {Zaf 
6 6 



4QQ 

-1.19568(1) + — ln[(Za)-^] 



(91) 



F(4P3/2) = ^ 



-lnA;o(3P) + (Za)2 



1 ^^7 

-0.63094(1) + ^ln[(Za)-2] 



(92) 



The results obtained for A^^ and Ag,! are in agreement with those previously known 0. 
The values of the Bethe logarithms P, 



In A;o(3P) = -0.03819023(1), In A;o(4P) = -0.04195489(1) 



(93) 



could be verified numerically with a 7 figure accuracy from our analytic expressions by 
numerical (Gaussian) integration. For Ag^i, we use the following general formula which may 
be extracted from the work by Erickson and Yennie ( [^], Eq. (4.4a) ibid., upon subtraction 
of the vacuum polarization contribution implicitly contained in the quoted equation): 



[1-5, 



ifi) 



(3-(/(/ + l))K) 



(2 / - 1) (2 /) (2 / + 1) (2 / + 2) (2 / + 3) 



+5i,i 



1-4 

601 



11. 

Y5 + i'^.v-i/2 



77 



240 60 n2 



71n2 + 3 (7 - Inn + ^(ra + 1)) 



(94) 



where 7 is Euler's constant, and ^ refers to the logarithmic derivative of the F function. For 
P states (Z = 1), this formula reduces to 



if3-i) + (i-i) 



15 



(95) 



and is in agreement with our results. From our results for F, we extract the following new 
values for the coefficients Aqq: 



^6,o(3Pi/2) = -1.14768(1), A6,o(3P3/2) = -0.59756(1), 



(96) 



and 



A,o(4Pi/2) = -1.19568(1), A6,o(4P3/2) = -0.63094(1). 



(97) 



The new results for ^0,0 are the main results of this work. They are in excellent agreement 
with data obtained from numerical calculations by one of the authors (PJM) and Y. K. Kim 
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0. Mohr and Kim calculate the F function defined in (|]) numerically for Z > 10, treating 
the binding field non-perturbatively. By extrapolating their numerical data to the region 
of small Z, we obtain the following estimates for the remainder function Gse,7 implicitly 
defined in Eq. (H). 



and 



G'sE,7(3Pi/2,^= 1) = 3.6±0.5, GsE,7(3P3/2,^ = 1) = 2.6±0.5, (98) 



G'sE,7(4Pi/2,^ = 1) = 3.9±0.5, G'sE,7(4P3/2,^ = 1) = 2.8±0.5. (99) 



The uncertainties in Gsej are used to estimate the theoretical uncertainty from the one- 
loop contribution. When modeling the numerical data, it must be taken into account that 



as noted by Karshenboim i coefficients vanish for P states. Values of Aq q and Gse,7 



for 2P states are given in |l|. We use the following implicit definition of the Lamb shift C: 



2 



E = mr [f{n,j) - 1] - — — ^ [/(n, j) - 1]^ + £ + Ems, (100) 

2[m + rriN) 

where E is the energy level of the two-body-system and f{n,j) is the dimensionless Dirac 
energy, m is the electron mass, rrir is the reduced mass of the system and mjy is the nuclear 
mass. It should be noted that we consider the hfs-fs mixing term as a contribution to the 
hyperfine structure. The small hfs-fs- mixing correction which is discussed in |^T|, mixes the 
F = 1 sublevels of the P1/2 and P3/2 states and shifts the center of the hyperfine levels. It 
should be taken into account when the fine structure is deduced from precision experiments. 

In order to calculate the Lamb shift, we include the Barker-Glover correction to hydrogen 
energy levels |2^, which we refer to as the [Za)^ recoil correction. We also include the {ZaY 
recoil correction calculated by E. Salpeter [^], and the results for recoil corrections of 
order (Za)^ rrir/mN- The results for recoil corrections of order (Za)^ rrir/misf have been 
confirmed by K. Pachucki (see ch. 5 of [|I|). We include contributions from the higher order 
two-loop correction of order (a/vr)^ {Zct)^ ln^(Za;)~^ corresponding to the -B6,2^coefficient 
P3| , and for three-loop corrections in lowest order. The theoretical error from the two-loop 
contribution (-Be.i and higher terms) is estimated as half the contribution from the recently 
calculated i?6,2^coefficient. The higher order contributions due to vacuum polarization of 
order a/vr (Za)^ can be obtained by analyzing the small distance behavior of the Dirac wave 
function, i.e. by evaluating the matrix element of the Uehling potential (see e.g. |2j 



a{Za) 

Vvp(r) ^ 



15 60 ^ 



(101) 



with P-state wavefunctions expanded in powers of Za. We obtain the results 



and 



3 ^2 _ 1 

Al^{nPy,) = --—^ (102) 



2 77^ — 1 

^60^(-P3/2) = -Y^^ (103) 
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for the leading term. We have also evaluated the contribution of the Uehling potential 
numerically without expansion in Za with the results 



Gu,7(3Pi/2, Z=\) = 0.0455, Gu,7(3P3/2, Z = V)= 0.0249, (104) 

and 

Gu,7(4Pi/2, Z = \) = 0.0480, Gu,7(4P3/2, Z = \) = 0.0262 (105) 

where the function G'u,7 is defined in analogy with Gsej- The contribution of the higher- 
order terms is negligible compared to the uncertainty in the higher-order self energy terms. 
The Wichmann-Kroll vacuum polarization contribution is expected to be of order iZa)'^ 
times the Uehling correction and is not included here. 

The above mentioned contributions to the Lamb shift are listed in the Tables |T| and |V 



for the states under investigation. It should be noted that the reduced mass dependence of 
the terms must be restored in low Z systems to obtain the correct value for the Lamb shift. 
Terms which are caused by the anomalous magnetic moment of the electron acquire a factor 
(m^/mg)^ (where niy is the reduced mass of the system and mg the mass of the electron), all 
other contributions to the Lamb shift acquire a factor im^jm^'^ . In addition, the argument 
of the logarithms ln[(Za)~^] must be replaced by ln[(me/mr Z a)~^]. The relevant formulae 
are also given in |^). 

It should also be noted that for two-loop and three-loop corrections in respective lowest 
order (a/vr)^ (^Q;)^ and (a/vr)^ (ZaY^ only the anomalous magnetic moment of the electron 
contributes to the Lamb shift for P states, because the Dirac form factor -Fi(g^) is infrared 
finite in two- and three-loop order. So it is only the contribution from the magnetic form 
factor -^2(5^ = 0) which persists. The calculation of the contribution to the Lamb shift is 
then straightforward. 

We obtain the following theoretical results for the Lamb shift of 3P and 4P states: 

/:(3Pi/2) = -3473.75(3) kHz, (106) 
/:(3P3/2) = 4037.75(3) kHz, (107) 
/:(4Pi/2) = -1401.52(1) kHz, (108) 

£(4P3/2) = 1767.30(1) kHz. (109) 

The theoretical values for the fine structure splitting, using the 1987 Cohen-Taylor value of 
a-^ = 137.0359895(61) |28|, are as follows: 

AEfs(3P) = 3250089.8(3) kHz (110) 

AEfs(4P) = 1371130.0(1) kHz (111) 

For 2P states, the theoretical value is A£;fe(2P) = 10969043(1) kHz 0. The uncertainty in 
the theoretical values for the fine structure splitting is given by the uncertainty in a. Any 
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determination of the fine structure beyond the quoted uncertainty would yield a value of a 



improved with respect to the 1987 value. Given the scattering of available data for a 



such a determination could be useful for checking the consistency of measurements coming 
from different fields of physics. 

The formula for the fine structure as a function of a for atomic hydrogen [Z = 1) is as 
follows: 



a 
+a' 

+a 



2n^ 



^2 y3 

2n^ 



+ a" 



2 vrra^ 



'a'^h' 7y 9y 3y 



xy 



32^3 16^4 4n5 An^(l + x) 



5/43) ^2y 



jl-3 ^3 



A. / X 1 -1 In^-l,,, 
AAe.o n + — — In y a 

lb 



-1 „ -2\ 



3 n 

2y^ , ^, , 31?/ 45v 7y 45?/ 15?/ 
AG{n) + + + ^rr^ - + ^ 



+ (112) 



256 n3 128 64 32^6 16 n^^ 
where the theoretical uncertainty in the difference of the Lamb shift of nP states is given by 
5Cti,{2) = 80 Hz, SCthiS) = 30 Hz, 5£th(4) = 10 Hz. (113) 
The mass ratios are 

x = me/mp and y = mr/rrie = 1/(1 + x). (114) 

AAe o(n) and AG{n) are defined as 

AAe,oin) = A,o(nP3/2) - A6,o(nPi/2), AG(n) = GsE,7(nP3/2) - GsE,7(r^Pi/2). (115) 

For practical purposes, the n-dependence of AGse,7(^) may be suppressed, because it is a 
very small contribution (in the 1 Hz range), and we may assume AGse,7('^) ~ —1.0. The 
two- and three-loop coefficients to the anomalous magnetic moment are given by [^| 



(e) 
02 



-0.328478965 and af^ 



1.18124156. 



(116) 



VI. CONCLUSIONS 

The analytic calculation of higher order binding corrections to the Lamb shift of excited 
P states has been described in this paper. We provide more accurate theoretical values of 
the Lamb shift for 3P and 4P states in hydrogenlike systems. 

We also give a formula for the fine structure as a function of a (for 2P, 3P and 4P states), 
which may be used to determine a from an improved measurement of the fine structure. 

With the possibility of substantial improvement in the precision of spectroscopic experi- 
ments (trapped atoms), a better determination of a from measurement of the fine structure 
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might be within reach in the near future. Such a determination of the fine structure con- 
stant a, from the effect on which its name is based, would complement other high precision 
determinations from solid state physics and the anomalous magnetic moment of the electron. 

We note that there are deviations of experimental data for excited nS-nP transitions 
from theory by more than one standard deviation but less than two standard deviations 
(see [0] and references therein). However, both theory of the Lamb shift and spectroscopic 
techniques have improved since the measurements were made, so one might expect a more 
precise comparison of theory and experiment in the future. The present uncertainty in the 
theory would in principle allow a determination of the fine structure constant with a relative 
uncertainty of less than 5 parts in 10^. However at this level of precision additional theoretical 
work might be needed to address questions such as asymmetries in the natural line shape. 
We only mention that for excited states, an experimental determination of the fine structure 
could be simplified by the slower decay (narrower line width) of the higher excited P states 
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TABLES 



contribution 


3-Pl/2 


3-P3/2 


Fnq 


-1.433010(1) + 248/405 In {e/{Zaf) + 


-1.433010(1) + 248/405 ln{e/{Zaf) 


FSy 


0.922653(1) - 20/81 ln(e/(Za)2) 


0.629717(1) -20/81 ln{e/{Za)^) 


FSH 


0.356318(1) - 73/324 In (e/(Za)2) 


0.333053(1) - 55/324 In (e/(Za)2) 


FSE 


0.0406519(1) + 1/36 In (e/(Za)2) 


0.013551(1) + 1/108 In {e/{Zaf) 


F54, 


-0.830340(1) + 40/81 ln{e/{Zaf) 


-0.236869(1) + 13/108 In {e/{Zaf) 


sum 


-0.94378(1) + 268/405 In {e/{Zaf) 


-0.69356(1) + 148/405 In (e/(Za)2) 


TABLE 1. 


Contributions of relative order {Za)'^ to the low-energy part Fl for the SP1/2 and 


3P3/2 states. 






contribution 


4Pl/2 


4P3/2 


F 

^ nq 


-1.512220(1) + 229/360 In {e/{Za)'^) 


-1.512220(1) + 229/360 ln{e/{Zaf) 


FSy 


0.966398(1) - 23/90 In {e/{Zaf) 


0.662154(1) - 23/90 In {e/{Zaf) 


FsH 


0.364541(1) -2891/11520 ln(e/(Za)2) 


0.342940(1) - 439/2304 In (e/(Za)2) 


FSE 


0.0335504(1) + 13/768 In (e/iZa)"^) 


0.012904(1) + 5/768 In (e/(Za)2) 


Fs., 


-0.850066(1) + 787/14-10 In (f/(Zn)2) 


-0.236345(1) + 53/288 In (r/(Zn)2) 


sum 


-0.99780(1) + 499/720 In {e/{Zaf) 


-0.73057(1) + 137/360 In (e/(Za)2) 



TABLE II. Contributions of relative order (Za)'^ to the low-energy part Fl for the 4Pi/2 and 
4-P3/2 states 
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J. •! J. • 

contribution 


7~> * 1 TT 

3-ri/2 in kilz 


7~> • 1 TT 

3-P3/2 m kHz 


one-loop sell-energy 


—3477.349(5) 


4046.413(5) 


two-loop self-energy 


7.705(23) 


-3.782(23) 


three-loop self-energy 


— L).L)o4 


0.032 


vacuum polarization 


— (J.lzz 


—0.027 


(^Zcx)^ recoil 


0.641 


-0.320 


(Za)^ recoil 


-4.705(13) 


-1.915(13) 


{Zaf recoil 


0.139 


0.139 


Sum for 3P 


-3473.75(3) 


4037.75(3) 


TABLE III. Contributions to the Lamb shift in kHz for the 3Pi/2 and 3P3/2 states. 


contribution 


4Pi/2 in kHz 


4P3/2 in kHz 


one-loop self-energy 


-1403.102(2) 


1770.887(2) 


two-loop self-energy 


3.252(10) 


-1.594(10) 


three-loop self-energy 


-0.027 


0.014 


vacuum polarization 


-0.054 


-0.012 


[Zoif recoil 


0.270 


-0.135 


[Zaf recoil 


-1.915(5) 


-1.915(5) 


{Zaf recoil 


0.061 


0.061 


Sum for 4P 


-1401.52(1) 


1767.30(1) 



TABLE IV. Contributions to the Lamb shift in kHz for the 4Pi/2 and 4P3/2 states. 
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FIGURES 



FIG. 1. Feynman contour for uj integration (one-loop self energy). Lines directly below and 
above the real axis denote branch cuts from the photon and electron propagator. Crosses denote 
poles originating from the discrete spectrum of the electron propagator. 
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FIG. 2. Mohr's contour for evaluating the one-loop self energy contribution to the Lamb shift. 

1 




FIG. 3. The cj-intcgration contour used by Pachucki and in the calculation presented in this 
paper. For the divergent terms in the high-energy part, we use the Wick-rotated contour given by 
the lines extending to e it zoo. For the naively convergent terms, we use the original contour Ch 
which extends to +00 rtid. ' 
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FIG. 4. Expansion of the bound electron self energy in powers of the binding field. 
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